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Abstract 



The paper addresses the question whether the dynamo remains of the solar 
type if it is coupled with the two-cell in radius meridional circulation pattern. 
The answer is positive for a wide class of dynamo models that take into account 
the subsurface rotational shear. We show that the Gleisberg-type cycles, repre- 
senting variations of 11-year on the centure time scale, can be generated due to a 
nonlinear resonance between the dynamo wave, and the large-scale magnetic field 
amplification in the middle of the convection zone as a result of the convergent 
meridional shearing flows. The conditions of such resonant interaction depends 
on the speed of the meridional flow and other details of the model. 

Subject headings: Dynamo — Magnetohydrodynamics (MHD) — Sun:dynamo 
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1. Introduction 

It is widely believed that the sunspot activity of the Sun is governed by the 
hydromagnetic turbulent dynamo operating deep in the solar convection zone. The 
meridional circulation, which is observed on the surface as a steady flow from the equator 
to the poles with a speed of 10-20 m/s, is often considered as an important ingredient of 



the solar dynamo (Dikpati & Charbonneau, 1999; Charbonneau, 2011). It is suggested 



that the return meridional flow in the deep convection zone can promote the drift of 
sunspot activity to the equator in the course of the 11-year solar cycle (so-called "butterfly 
diagram"), as well as the meridional drift of the poloidal magnetic field on the solar surface 
from mid-latitudes to the poles. This idea is extensively used in the Babcock-Leighton type 
solar dynamo models. Despite, the mean-field dynamo models may reproduce the main 



properties the solar cycles without the meridional circulation (e.g., Pipin & Kosovichev 



2011b, Pipin 2013), the circulation affects details of the model. In particular if the dynamo 



effect is distributed in the bulk of the convection zone, the circulation properties affect the 
magnetic field distribution with depth. For example, the meridional circulation pattern, 
which is derived from a mean field model of the solar differential rotation has the latitudinal 



flow concentrated to the boundaries of the convection zone (Kitchatinov, 2011). Such 



circulation swaps the toroidal magnetic fields to the bottom of the convection zone, and 



the poloidal magnetic fields get concentrated to the pole (Pipin & Kosovichev, 2011a). The 



helioseismology inversions suggest that the global circulation may consist of two cells in the 



radial direction (see, Mitra-Kraev & Thompson 2007, Zhao et al. 2012). In the northern 



hemisphere, the flow direction is from the equator to the poles at the top of the convection 
zone, and the flow direction is reversed in the bottom cells. Such double-cell circulation is 



also obtained by numerical simulations (Miesch et al. 2006, 2011; Guerrero et al. , 2013). 



The influence of multi-cell circulation on the flux-transport dynamo models was studied by 



Bonanno et al. (2006); Jouve & Brun (2007). They found that in case of the double-cell 
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circulation pattern, the evolution of the large-scale toroidal and poloidal magnetic field of 
the Sun is qualitatively different from the observations. 

In this Letter we examine effects of the double-cell meridional circulation pattern (see, 
FigHJ) on a mean-field dynamo model, in which the magnetic field generation is distributed 
in the bulk of the convection zone, and the migration of the magnetic field on the surface is 
controlled by the subsurface rotational shear layer action. This type of the subsurface-shear 



shaped dynamo was suggested by Brandenburg (2005) and was extensively studied in our 



recent papers (e.g., Pipin &: Kosovichev||2011c ). In particular, it was found that the dynamo 



model can reproduce the known statistical relation of the solar cycle (Pipin et al., 2012). 



Here, we investigate effects of the double-cell meridional circulation pattern on this model, 
and show that the resulted evolution of the large-scale magnetic field is in qualitative 
agreement with observations. 



2. Dynamo equations 

We consider a large-scale axisymmetric magnetic field, B = e^B + V x A ^ , where 
B(r,9,t) is the azimuthal component of the magnetic field, A(r, 9,t) is proportional to the 
azimuthal component of the vector potential, r is the radial coordinate, and 9 is the polar 
angle. The mean flow is given by velocity vector V = e r U r + eglJg + e^r sinOQ, where 
Q(r,9) is the angular velocity of the solar differential rotation, and U r (r,9) and Ug(r,6) 
represent velocity components of the meridional circulation. The mean magnetic field is 
governed by the dynamo equations: 

¥ =Vx(VxB + £), (1) 



where mean electromotive force £ is given by [Pipm] fl2008| ) (hereafter P08) : 



Si = fa + 7 ^ A) ) B 3 - fa + 77$) V s B h . (2) 
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( H) 

The tensor a^, representing the so-called a-effect, includes hydrodynamic (a\j ) and 
magnetic (a£^) helicity contributions: 



(H) 



C a ip a sin 2 9a\j + a\j , (3) 

h fa (/# (e • A«) + A ( ? (e • A<»>) ) } + (4) 

+ eiei {3^(/W(e.AW)+/W(e.AM))} + 

3r? T { (e^W + e . A W) ;W + + e . A W) /g W} ; 

where A^ = Vlogp quantifies the density stratification, A^ = ^Vlog [VtJ quantifies 
the stratification of turbulent diffusivity 77^? , and e = ft/ \ ft | is a unit vector along the axis 
of rotation. The a-quenching function ip a = —3/40^ depends on (3 = \B\ / \Jpfi u 2 , 
and 0g is given in P08. 

The turbulent pumping, %f \ depends on the mean density and turbulent diffusivity 



stratification, and also on the Coriolis number Q* = 2t c Q , where r c is a typical convective 
turnover time, and Qq is the angular velocity. Following the results of P08, 7^ is expressed 
as follows: 

7^ = 3 VT {ft ) Al p) + fi a) (e-\^)e n }e m3 -3 VT fi a) e 1 e mm e n A^, (5) 
- 3 VT (e - 1) {/ 2 (a) A^ + f[ a) (e ■ A«) e n ) e inj . 
The turbulent diffusivity is anisotropic due to the Coriolis force, and is given by: 

Vijk = 3r] T I (2f[ a) - / 2 (d) j e ijk - 2fi° :) e i e n e njk j . (6) 

We also include the nonlinear effects of magnetic field generation induced by the large-scale 
current and the global rotation, which are usually called as the Q x J-effect or the § dynamo 



effect (Radler, 1969). Their importance is supported by the numerical simulations (Kapyla 



et al. 


2008; 


Schrinner, 


2011 



2011). We use the equation for r)fX which was suggested in P08 



(also, see, Rogachevskii & Kleeorin, 2004): 



Vijk - 6 Vt^&Ja e j i ^7 d ik + ^2 



~(w) BiB k 



B 



(7) 
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where, C$ measures the strength of the Q x J effect, if 27 (P) are normalized versions 
of the magnetic quenching functions ip 2 W j given in P08. They are defined as follows, 



^27 (P) = fv 9 ^ (ft)- The functions /j^^m Eqs(3 5 6 7) depend on the Coriolis number. 



They can be found in P08, as well. 



In our models we use the solar convection zone model computed by Stix (2002). 
The mixing-length is defined as £ = ckmlt 

|A^)| , where = Vlog p is the inverse 
pressure scale height, and «mlt = 2. The turbulent diffusivity is parameterized in the form, 



VT = Crff, where ^ 



(0) 



U ' 2 T r 



is a profile of the mixing-length turbulent diffusivity, £ is 



3/™ (r) 

a typical correlation length of the turbulence, C v is a constant to control the efficiency of 
large-scale magnetic field dragging by the turbulent flow. Also, we modify the mixing-length 
turbulent diffusivity by factor / OT (r) = 1 + exp (50 (r ov — r)), where r ov = O.725i? to 
model saturation of the turbulent parameters near the bottom of the convection zone, as 



Ossendrijver et al. 


2001, 


2002; 


Kapyla 



et al. 2008). The results do not change much if we scale the turbulent diffusivity gradient 
\( u ) — Q v Vlog \ VtJ w ffh a factor C v < 0.5. For values of greater C v our model leads to a 
steady non-oscillating dynamo which is concentrated to the bottom of the convection zone. 



Note, that the models by Bonanno et al. (2002); Bonanno (2012) use C v = 0. We introduce 



the additional parameters C v and function f v{r) to model the distribution of the helicity 
coefficient a closer to the results of the numerical simulations. 

The bottom of the integration domain is at r& = O.715_R , and the top of the integration 



domain is r e = 0.99i? Q . Pipin & Kosovichev (2011a) showed that the solar-type dynamo 
solution can be obtained for C a /C$ > 2. In that paper we found that the approximate 
threshold for magnetic field generation is C a ~ 0.03 for a diffusivity scaling factor C v = 0.05 
to match the 22-year periodicity. Figure [T] shows the radial profiles for the a-effect 
coefficients, the isotropic and anisotropic parts of the turbulent diffusivity and the Q x J 
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effect. They are in the qualitative agreement with the simulations results by |Ossendrijver 



et al. 


(2001 


) and 


Kapylii et al. 


(2009) 



The meridional flow is modeled in the form of four stationary circulation cells with two 

cells along the radius in each hemisphere. The pattern is modeled by the stream function 

U = ^Vx (e^), where we assume the following analytic form to describe the radial 
P 

and latitudinal dependence of 



2c m / 2 2\ ( 1 

-K \X 



1.5 



sin 



7T (X — Xi 



lx e 



X 



' dP dP 

x (tanh (20fc) (x 2 - xj - l) + l) ( + m '' 



36 89 

where, Uq is the amplitude of the flow, x = r/R is a distance in the units of the solar radius, 
Xm — r m/R is the stagnation point of the upper cell, X{ = Ti/R is the inner boundary 
of the integration domain, x e = r e /R is the upper boundary of the integration domain; 
parameter m controls the number of cells in latitude, and parameter k is to control the 
strength of the flow at the bottom of the convection zone; c m is the constant to normalize 
the maximum of the flow amplitude to 1. In the model, we use x m = 0.89, k = and 
m = —1/3, that gives c m = 1/14 for Xi = 0.715 and x e = 0.99. The given parameters 
result in the circulation pattern shown in Figure lc. This pattern resembles qualitatively 



the results of the helioseismology inversions (Mitra-Kraev & Thompson, 2007, also Zhao et 



al 2012). Note, that for x m = and keeping the others parameters the same we would get 
a single-cell (per hemisphere) flow pattern which is close to the mean-field models results 

(see, 



Kitchatinov 2011). If we additionally put k — 1, then the amplitude of the flow near 
the bottom of the convection zone is reduced strongly. For m > 1/3, secondary latitudinal 
cells become pronounced at the high latitude . 



The contribution of small-scale magnetic helicity \ = a • b (a and b are the fluctuating 
parts of magnetic field vector-potential and magnetic field vector) to the a-effect is defined 



-8- 



as follows (P08) 



(a)r XTc n Act) XJc 

fi p£ 2 fioP^ 2 



= 2ft a >5. 



(9) 



The nonlinear feedback of the large-scale magnetic field to the a-effect is described by a 
dynamical quenching due to the constraint of the total magnetic helicity conservation. The 



latter can be written as follows (Hubbard & Brandenburg, 2012, Pipin, 2013): 



d 
dt 



J {x + A-B}dV = -nj {B- J + b-j}dV 



(10) 



where integration is performed over the volume that comprises the dynamo region, J-* is 
the diffusive flux of the total helicity which results from the turbulent motions, rj is the 



molecular diffusivity. The differential equation that corresponds to the Eq.(10) is: 



dx 
dt 



d(A-B) 

dt 



X 



R T 



-7]BJ-V-T X - (U- V) (x + A-B). 



'IV, 



where J- x = — ?7 X V (x + A ■ B). We assume R m = 10 6 and rj x = O.lrjT- To solve Eq.([TT]) 
the large-scale vector-potential A is required. For axisymmetric large-scale magnetic fields 
the vector-potential can be represented as a sum of the azimuthal and poloidal components 



(Krause & Radler, 1980): 



A = e/T + rP 



r sin 9 



A + re r P. 



12) 



The toroidal part of the vector potential is governed by the dynamo equations (Eq|TJ). 
The poloidal part of the vector potential can be restored from the toroidal magnetic 
field component using equation V x (rP) = e^B. The reconstruction is simple for 
the pseudo-spectral numerical schemes which are based on the Legendre polynomial 
decomposition for the latitude profile of the large-scale toroidal field. We matched the 
dynamo solution to the potential field model outside the convection zone, and assume the 
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perfect conductivity condition at the bottom boundary. For the magnetic helicity we employ 
X = at the bottom of the convection zone and we assume that the the radial derivative 
of the total helicity is zero at the top. The initial magnetic field is assumed dipolar with a 
small quadrupolar component to check the parity preference when the solution reach the 
steady state. 

3. Results and discussion 

Figure [2] shows the snapshots of the magnetic field evolution in the North segment of 
the solar convection zone for the models without and with the meridional circulation. We 
observe drift of the dynamo waves related to the evolution of the large-scale toroidal and 
poloidal fields, towards the equator and towards the pole, respectively. It is found that 
circulation accelerates the equator-ward drift of the toroidal magnetic field near the surface 
and the poleward drift of the toroidal magnetic field the near bottom of the convection 
zone. The strength of the radial magnetic field at the poles is increased as well. It is 
also found that the double-cell circulation reduce the toroidal magnetic field strength 
inside the convection zone. This effect is generally expected for the circulation-dominated 
flux-transport dynamo models, for which an increase of the circulation speed increases the 
dynamo generation threshold in terms of the a-effect coefficient. This, in turn, decreases of 
the period of the dynamo cycle. 

The time-latitude "butterfly" diagrams of the toroidal magnetic field evolution in the 
upper part of the solar convection zone and the similar diagrams for the radial magnetic 
field at the surface are shown in Figure [3j In all cases (with and without circulation) there 
is the qualitative agreement with observations. We found that the maximum of the toroidal 
magnetic field migrates closer to the equator for the models with Uq = 10 and Uq = 15 m/s. 
However, in these cases the butterfly diagram wings are wider in latitude than in the case 
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9 = 45°; panel (b), the radial profiles of the background turbulent diffusivity C^rj^ , the 
isotropic, rj^\ and anisotropic, vf A \ parts of the magnetic diffusivity coefficient, and the 



Q x J term (Eq. un; (c), velocity vectors of the double-cell meridional circulation pattern . 




Fig. 2. — Snapshots of the magnetic field evolution inside the North segment of the convection 
zone for the dynamo model: (a) shows the model without meridional circulation; (b), the 
model with the meridional circulation with the maximum speed Uq = 10 m/s. The field lines 
show of the poloidal component of the mean magnetic field, and the toroidal magnetic field 
(varies ±lkG) is shown by gray scale background images. 
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without circulation. Also it is found that the circulation reduces the latitudinal width of the 
polar branch for the radial magnetic field evolution and also reduces the overlap between 
the cycles. We draw these diagrams only for one hemisphere because in all these cases the 
antisymmetric mode (dipolar) is dominant. 

However, the changes of the dynamo properties with increase of Uq from to 15 m/s are 
not monotonic. We found a resonance-like phenomenon for the model with the flow speed 
Uq — 5 m/s. For this case, the evolution of the magnetic fields shows long-term periodic 
variations, which correspond to a torus-like attractor on the dynamo phase diagram (see, 



e.g., Knobloch et al. 1998). In this case, the short-term (11-year) cycles are modulated by a 
grand cycle with a period of 80-90 years, corresponding to the so-called "Gleisberg" cycle. 
Figure p] shows some results for the model with Uq = 5 m/s. We found that the Gleisberg 



cycle appears irrespective of the magnetic helicity feedback (Eq. 11). The only difference 
is that for the model without the magnetic helicity feedback, the long-term variations keep 
the dipole-like symmetry relative to the equator. In the model with the magnetic helicity 
feedback, the maximum of the magnetic activity periodically drifts from one hemisphere 
to another on the long time scale. Figure p] shows the simulated sunspot number (defined 



by Pipin et al. 2012) and the polar flux variations. The polar flux is defined as the 



total B r above ±60° latitude. The model shows that the polar flux reversal time can be 



different in the two hemispheres. The similar effect was reported earlier by Pipin ( 1999 ) for 
the nonlinear model that included a nonlinear feedback of magnetic field on the angular 
momentum balance in the solar convection zone. Figure [4](c) shows that the average polar 
flux in the northern hemisphere is larger than in the southern one. This is because of the 
asymmetric initial conditions and due to the parity conservation in the model. 

We can interpret the obtained long-term variations as a nonlinear resonance between 
the dynamo wave propagation and the increase of the magnetic field in the middle of the 
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Uo = m/s 





TIME, [YR] 

Fig. 3. — (a) the time - latitude variations ("butterfly diagrams") for the model without 
meridional circulation; (b) the same for the model with the meridional circulation with the 
maximum speed Uq = 10 m/s; (c) the same as (b) for the Uq = 15 m/s. The toroidal field 
near the surface, at r = O.95_R , is shown by contours (plotted for ±100G range) and the 
radial magnetic field at the surface is shown background images. 



convection zone due to the convergence of meridional flows from the top and bottom of 
the convection zone. The long-term modulation (Gleisberg cycle) disappears if we discard 
the nonlinear part of the 5 effect (ftx J term in Eq|7]). Note that in Eq.Q the function 
ip^ (ft) is monotonically quenched to zero when the magnetic field equipartition parameter 
ft is growing. Function (jffi (ft) is non-monotonic, and has a maximum at ft « 1. Moreover, 
the structure of the second term in Eq. ^ is the same as for the helicity coefficient a effect 
because the associated electromotive force is proportional to l^rC 'g fjjf' '^2"' B '% ( e ' ^) logB 2 



see details in P08 and also in Rogachevskii & Kleeorin 2004). Our interpretation for the 



origin of the long-term modulation is supported by the fact that for the given dynamo 
number C a the long-term oscillations are excited only within a narrow range of the 
circulation amplitude, e.g., for C a = 0.025, the resonance circulation speed is f/o ~ 5 m/s. 
The long-term periodicity disappears for Uq ~ 3 m/s and Uq = 7 m/s. Moreover, for the 
Uq = 10 m/s the similar long-term periodicity is excited for C a = 0.035. We guess that 
the similar resonance could exist for the mean-field dynamo with double cell meridional 



circulation and the a effect due to buoyancy instability (see, e.g., Chatterjee et al. , 2011), 
because the latter is proportional to V-B 2 as the nonlinear part of the 5 effect. 

For the excitation of the North-South asymmetric grand cycle the effect of the magnetic 
helicity is important. The growth rate of the quadrupolar modes is smaller than for the 
dipolar ones. This leads to the quenching of the asymmetric perturbations. The nonlinear 
feedback of the magnetic helicity on the magnetic field generation results in the time-delay 
of the a-effect quenching. When the typical time-delay become comparable with the growth 
time of the grand cycle the asymmetric perturbations leads to excitation of the long-term 
asymmetric cycles. The nonlinear time-delay feedback on the magnetic field generation 



is usually considered as important mechanism of the grand cycles excitation (Yoshimura 



1978, Kichatinov et al. 1994, Covas et al., 1998, Knobloch et al., 1998 Pipin, 1999). 
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Fig. 4. — Panel (a) toroidal (contours) and radial (background images) fields for the speed 
Uq = 5 m/s; panel (b) shows the simulated sunspot number (SN) as a function of time (see, 



Pipin et al. 2012); (c) the corresponding polar flux in the northern (solid) and southern 



(dashed) hemispheres. 
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We investigated the stability of the obtained solutions with respect to changes of the 
circulation pattern and the other parameters of the model. In addition to the previous 
discussion, we found that the long-term modulation disappears when we neglect the 
convective overshooting effect, i.e., when we put f ov (r) = 1. The results similar shown 
in Figures 2(b) and 3(b,c), are obtained for the case of pure a a 2 Q dynamo (without the 
Six J effect). However, in this case, the cycle overlap in the time-latitude diagrams is much 
stronger than in Figures 3(b,c). The similar time-latitude butterfly-diagrams are obtained 
also for the case when the meridional circulation forms only one cell per hemisphere with 
the stagnation point radius larger than 0.85-R. This is because of the subsurface rotational 
shear layer included in the model. For the deeper stagnation point the dynamo wave is 



dominated by the toroidal magnetic field near the bottom of the convection zone (Pipin 



& Kosovichev, 2011a). Also, we checked the effect of the meridional circulation in the 



simplified model which was presented by Pipin & Kosovichev (2011c). In that model the 



ct-effect was confined in the low-latitude ±30° zone. The results are generally similar to 
what is shown in Figures 2(b) and 3(b,c). 

The main findings can be summarized as follows. The mean-field dynamo model that 
includes the subsurface rotational shear layer and the double-cell (in radius) meridional 
circulation, indicated by the recent helioseismology results, can reproduce the solar magnetic 
cycles in the form of the time-latitude "butterfly" diagrams. The double-cell circulation 
affects the distribution of the magnetic field with radius in the convection zone, increasing 
the field concentration to the convection zone boundaries, and to the layer of the two cells 
convergence in the middle of the convection zone. The latter effect can lead to a result in a 
nonlinear resonance between the dynamo wave accompanied by large-scale magnetic field 
amplification in the middle of the convection zone as the result of the convergent flows. 
The resonance conditions depend on the flow speed and other details of the model. This 
resonance leads to a long-term modulation of the dynamo cycles, resembling the Gleisberg 
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100-year solar cycles. In this case the solution can be North-South asymmetric due to 
unstable quadrupole component. 
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